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Abstract 

A new Monte Carlo implementation of Djordjevic's dynamical scattering generalization of the DGLV radiative 
energy loss opacity series is used with a hybrid interpolation scheme to compute both light and heavy quark jet 
quenching up to third order in opacity. The enhancement of the ratio of bottom to charm quark energy loss due to 
perturbative long range color magnetic effects in nonuniform Bjorken expanding geometries is found to reduce the 
significance of the heavy quark jet puzzle posed by the observed near equality (within sizeable errors) of pion and 
nonphotonic electron nuclear modification at RHIC. Jet Flavor Spectroscopy discussed below will be a powerful tool 
to differentiate competing dynamical models of the QGP produced in ultra-relativistic nuclear collisions. 



1. Introduction 

We report results of a recent detailed Monte Carlo studyfT] of dynamical magnetic scattering effects on the radia- 
tive energy loss of light and heavy quark jets in Au+Au collisions at ^fs - 200 AGeV. The aim of this work is to 
advance the theory of jet tomography for quantitative applications to heavy ion collisions in accord with the objectives 
of the DOE topical JET collaboration projectQ. 

Since 2005 the heavy quark jet puzzle|l3l has remainedpT] one of the outstanding unsolved problems within the 
perturbative QCD multiple collision theory framework. Predictions of heavy quark quenching based on the DGLV 
opacity series model|5| and other pQCD models of radiative energy loss were found to significantly over-predict the 
midrapidity nonphotonic electron yield in the pr ~ 5 GeV range due to the relatively small energy loss of heavy 
bottom quark jets. Even after inclusion in WHDGlSl of additional elastic energy loss and fluctuations of path lengths 
effects (as required to maintain consistency of the theory with the observed light parton/pion nuclear modification 
factor R'^^ ~ 0.2), the bottom quark contribution to the nonphotonic electron spectrum was found to be too weakly 
quenched to explain the PHENIX and STAR data|3|. This discrepancy has led to novel gravity dual holographic 
models of jet energy loss completely outside the pQCD framework (see |4 7 8 1 and refs. therein). 

In this work we focus only on recent refinements of perturbative Hard Thermal Loop (HTL) QCD dynamical 
theory that have not been quantitatively studied up to now. In particular, we follow the proposal of (MD) M.Djordjevic 
||9l that long range dynamical color magnetic scattering in the HTL framework may significantly enhance the bottom 
quark radiative energy loss and thereby possibly help solve the heavy quark puzzle without having to invoke new 
non-perturbative dynamical assumptions. We have therefore developed a Monte Carlo hybrid interpolation version of 
the DGLV-BFW-MCLO code|10| to evaluate numerically the sensitivity of the ratio of charm and bottom quark jet 
quenching to perturbative color magnetic effects up to moderate order in opacity. 
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The hybrid DGLV-BFW-MC approach, used here, involves replacing the GW lfTTl (static Debye screened) mo- 
mentum distributions in the GLV(eq.l 13)/DGLV(eq.l7) multiple collision kernel 

with a path dependent effective magnetic enhanced transverse distributions 

_2 ^ i^e(zfN(r,„) 1 

^(Z,q,r„,) ^ ((q2+^2(z)))((q2 + 4^2(z))) ' 

where < r„, = fXm/f^e < 1 is the ratio of the color electric Debye and the assumed longer color magnetic screening 
lengths. The normalization factor Nir) = (1 - r^)/Log[l/P'] reduces to unity when r = 1 but has a weak logarithmic 
zero for r — > 0. This zero is however cancelled by the weak divergence of the effective elastic cross section o-(z, r,„) - 
9na^ /(2^^(z)N(rm)) as shown in detail by Djordjevic|9| for = 1. However, in higher opacity orders the unitarity 
corrections diverge in the (unphysical) r,„ - limit and require more care. Fortunately, in the finite size quark gluon 
plasma produced in A+A, r,„ is bounded from below by l/ifigRA) due to confinement of color outside the finite size 
plasma. We consider here finite r„, > 1/3 to explore the sensitivity of the bottom to charm jet ratio to enhanced soft 
momentum transfer due to perturbative magnetic field fluctuations of range up to 3/gT. 



2. Static Brick Problem 

At first order in opacity, the MD (r,„ = 0) modelf9l of the induced gluon number radiated per light cone momentum 
fraction from a massive quark jet produced at position x with energy E in azimuthal direction (p is given by 

^(x0) = ^ r^^^'^^'g 1 2(k+q) 
dx+ ' rr J 7T n x+ q^(q^+fp-) {k+q)^+x 

(2) 

where x+ is the fraction of plus momentum carried away by the radiated gluon, fi = gT{x, <p, t) is the local path 
dependent Debye screening mass and A^^,^ - 3asT(x,<p,T) is the inverse local dynamical mean free path. Here 
X = M^xl + OTj(l - x+) controls the "dead cone" effect due to the finite jet quark mass M as well as the local effective 
gluon thermal mass nig - ^^"^'^^ . For the idealized static "brick" problem Pqcp - 0{L - t). For r,„ > we use our 
interpolation distribution eq.(l). We include effects of Poisson fluctuations of the radiative energy loss due to gluon 
number fluctuations to compute P(e), the probability distribution of radiating a fraction of energy e. 

The left panel of Fig.l shows that the ratio of charm to bottom energy loss at first order in opacity is reduced by 
approximately 20% in the dynamical MD|9| compared to the static GW| 1 1 1 case, but the ratio for L > 2 fm remains 
well above unity in both limits. In the right panel the ratio of light (up) quark energy loss to bottom quark energy 
is shown vs path length 2 < L < 4 fm for different r„, and different orders in opacity. It is found that charm and up 
quarks have essentially identical energy loss. What is clear is that neither long range magnetic scattering nor higher 
orders are sufficient to bring the ratio near unity for L> 2 fm. 

Despite the fact that the absolute value of hE{Q) is almost doubled by switching from the pure static r,„ = 1 to the 
dynamical r,„ - model in the uniform brick geometry, as seen in Fig.2, the c/b and u/b ratios in Fig.l are surprising 
insensitive to the magnetic screening length. 



/ (k+q) k W U^+qf+X W , . 



3. Bjorken expansion 

To test the effect of magnetic enhanced scattering rates in a more realistic diffuse transverse geometry with lon- 
gitudinal Bjorken expansion cooling we assumed Wood-Saxon participant density profile and tq/t proper time de- 
pendence for pqgp{x.,t) in eq.2. We terminate the time integral when the Debye screening mass drops below Aqcd- 
Expansion reduces the effective path length and this helps to reduce the ratio of charm to bottom energy loss as seen 
from Fig.l. We average initial jet production points according to the Glauber Taa{x, b) profile and average over the 
jet angles. 
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Figure 1: Left: The ratio of charm to bottom jet radiative energy loss versus path length in a uniform static QGP at T = 250MeV 
is shown for E = 10 and 20 GeVfor the static G WUTl distribution (r,,, = 1 in eq.l) compared to the dynamic MDl^9j distribution 
(I'm = 0). Right: The ratio of light up quark (thermal mass) to heavy bottom quark energy loss vs L for E = 20 GeV jets and different 
r^ = 1,1/2, 1/3, 0. Solid lines correspond to N = I while dashed lines correspond to hybrid MC results with N = 1+2 + 3. The 
ideal r„ = case is the MD model[^. 



Analogously to the brick setup, for each of the jets we compute the probability distribution P(x, 0; e) and then 
average the results over the initial jet production profile in order to obtain the radiative energy loss AE/E and the 
partonic Raa- 

/ deeP(x,4>;€) Raa^ j ^dxpj„(x) j de(l - e)" P(x, cf>- e) . (3) 

The approximation in ([3]l is strictly valid only if the initial Au+Au spectrum is equal to the p+p spectrum and 
proportional to a constant index power law n, here set to 4. However, it provides a simple estimate that allows us to 
compare the MD and DGLV models in a static uniform cylindric geometry or a Bjorken expanding one. The results 
are shown in Fig 2, where we plot AE{M)/E and Raa(M) as a function of the mass of the quark for E - 2QGeV. 

The left panel of Fig. 2 shows that the dynamical MD scattering increases the energy loss fraction in a static plasma 
by about a factor of two. However, after including the dilution effects due to the expanding geometry, AE{M)/E 
decreases back toward the uniform static geometry, static GW value. Also quite remarkably, there is an approximate 
quark flavor independence from < M <2 GeV. Thus, identified charm quark jets can serve as calibrated light quark 
jet probes and hence can serve to uncover the gluon jet contribution to the pion Raa- Regarding the bottom quark, 
we find a modest increase of the energy loss with respect to GW and static geometries, but not sufficient to enhance 
AE(b)/AE(c) beyond about 0.5. 

On the right panel, /?aa(G) is shown after adding the mean elastic energy loss contribution from |6 | for a - 
0.3. The figure provides a convenient and intuitive "Jet Flavor Spectroscopy" representation of theory as well as 
experiment. It (1) emphasizes thatpQCD tomography with light (u,d,s) and charm quark jets is essentially equivalent, 
and thus tagged charm tomography fixes the light quark component of the high pT pion quenching Raa(^) ~ fgRAAig)+ 
fudsRAA(c) and (2) shows the mass splitting of Raa between b and c quark jets which depends on both the detailed 
pQCD dynamics as well the diffuse expanding QGP geometry. 

The shadowed regions represent the current rather large error bars in the pion and nonphotonic electron data. Once 
we take into account such large error bands, the inclusion of color magnetic effects in diffuse expanding geometries 
does help to reduce the significance of the heavy quark puzzle as posed in WHDG|6|. In this representation, the jet 
flavor spectroscopy hierarchy predicted in the pQCD framework including dynamical scattering as well as dynamical 
geometry 

RAAig) < RAAin ~{g + m)/2) < Raa(c) < RAA(e- ~ (c + b)/2)) < RAA(b) < RAAiy) 



may actually not be incompatible with present light and heavy jet tomographic data. 
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Figure 2: Illustration of Jet Flavor Spectroscopy. Left: Comparison between GW UIV and MDf^ average radiative energy loss 
AE(M)/E as a function of M, for uniform cylindric and Bjorken expanding geometries. Right: radiative plus elastic plus gluon 
fluctuation R/y^{M) vs parton flavor mass is shown for Bjorken expanding diffuse Glauber geometry including MD magnetic en- 
hancement and Os = 0.3. The empty circles represent the partonic contribution of the up,charm,bottom quarks. The filled circles 
are estimates of the pion and electron i?^^ assuming equal contribution of light quarks/gluons and charm/bottom jets respectively. 
The broad rectangular shaded box represents the current PHENIX and STAR data. Elimination of the wide horizontal uncertainty 
bar in RAA{e^) requires future c and b jet flavor tagged measurements. 

4. Conclusions 

Our main conclusion is that flavor tagged charm and bottom quark jet data is essential to confirm or reject the 
predicted pQCD nuclear modification flavor hierarchy that can be conveniently revealed using the "Jet Flavor Spec- 
troscopy" diagram Fig. 2b. The mass splitting between c and b jets is a particularly robust observable to differentiate 
pQCD and holographic gravity dual modelling of jet-plasma interactions, as emphasized in fT\. Our preliminary re- 
sults with hybrid Monte Carlo and 1/3 < r,„ = A'm/A'c < 1 up to third order in opacity indicate that magnetic effects 
help drive R^a down toward the present RHIC data and also reduce the b/c mass splitting. However, this hybrid pQCD 
approach still predicts that RAA(e^) / Raa(^) ~ 2 due to the highly quenched gluon jet component of the n yield ||6l. 
Once Raa(c) is known, the gluon Raa{<i) also becomes experimentally accessible since c ~ m in this framework. As a 
final point, we assumed no initial state saturation or shadowing eff'ects. Those have to be determined from p + A or 
direct photon Raa(7) to allow deconvolution of initial from final state quenching effects. 
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